Effects of nuclear absorption on the A/p ratio in relativistic heavy ion collisions 
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An enhanced A/p ratio in heavy-ion relative to p + p collisions has been proposed as one of 
the signatures for the Quark-Gluon Plasma (QGP) formation. A significantly large (A + £° + 
l.lE~)/p ratio of 3.5 has been observed in the mid-rapidity and low transverse momentum region 



in central Au+Au collisions at the nucleon-nucleon center-of-mass energy of 



4.9 GeV at 



the Alternating Gradient Synchrotron (AGS). This is an order of magnitude larger than the values 
in peripheral Au+Au collisions and p + p collisions at the corresponding energy. By using the 
Ultra-relativistic Quantum Molecular Dynamics (UrQMD) transport model, we demonstrate that 
the observed large ratio can be explained by strong absorption of p's (~99.9%) and A's (~99%) in 
dense nuclear matter created in central collisions. We find within the model that the initial A/p 
ratio, mainly from string fragmentation, does not depend on the collision centrality, and is consistent 
with that observed in p + p collisions. This suggests that the observed large (A + E° + l.l£~)/p ratio 
at the AGS does not necessarily imply the formation of the QGP. We further study the excitation 
function of the ratio in UrQMD, which may help in the search and study of the QGP. 
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PACS numbers: 25.75.-q, 25.75.Dw 

Strangeness production has been extensively studied 
in heavy-ion collisions because enhanced strangeness pro- 
duction may signal the formation of Quark-Gluon Plasma 
(QGP) This is due to the fact that the strangeness 

(ss) production threshold is significantly lower in a QGP 
than in a hadronic gas in which a ss pair has to be pro- 
duced by a pair of strange hadrons. Strangeness enhance- 
ment is often studied by the charged kaon production 
rate and the kaon to pion yield ratio (K/ir). All exper- 
imental results showed an unambiguous enhancement in 
kaon production rate and K/tt ratio in heavy-ion colli- 
sions with respect to elementary p+p collisions ■ 
However, the enhancement results can be also explained 
by particle rescattering as implemented in many hadronic 
transport models [9(. 

At AGS energies, the collision zone created in cen- 
tral heavy-ion collisions is baryon dense [Tfj| - [r3 |. In a 
QGP with high baryon density, production of light anti- 
quarks (u, d) should be suppressed, hence the ratio of 
anti-lambda to antiproton yields (A/p) should exhibit 
a larger value of enhancement than does the K/n ra- 
tio 0, EH, EH- This is because A/p gets enhancement 
not only from enhanced strange antiquark production but 
also from the suppressed it production. 

The E864 Collaboration at the AGS has deduced 
a (A + £° + l.l£~)/p ratio in Au+Au collisions at 
the nucleon-nucleon center-of-mass (c.m.s.) energy of 
y/s NN = 4.9 GeV. They deduced the ratio at mid-rapidity 
and almost zero transverse momentum (pr) by contribut- 
ing the discrepancy between th eir p measurement [15l . [16[ 
and that from AGS/E878 [13, 13 entirely to the differ- 
ent acceptances of the two experiments for p's from A 
decays [l! QJ]. The deduced (A + £° + \.\Tr)/p value 
has a strong dependence on the collision centrality. In 
peripheral collisions, it is consistent with p + p results 
(~ 0.2) at similar energies 0, |2(j- In most 10% central 



collisions, the (A + £° + l.l£~)/p ratio reaches a most 
probable value of 3.5, which is over an order of magni- 
tude larger than those in peripheral collisions and in p+p 
collisions. The E917 Collaboration has made direct mea- 
surements of p and A + E° yields at mid-rapidity and in- 
tegrated over px in central and peripheral collisions [2l| . 
The ratios of (A + S°)/p, given the large error bar, are 
consistent with E864. These results are intriguing be- 
cause they may point to possible QGP formation at the 
AGS. 

There are at least two physical origins for the large 
(A + S° + 1.1£~)/|> ratio: (1) an enhanced A/p (and/or 
E°' _ /p) ratio at the initial production stage of the an- 
tibaryons, and (2) a strong absorption of p's and a less 
strong absorption of A's and E's in nuclear matter pro- 
duced in heavy- ion collisions. An enhanced A/p ratio at 
the initial production stage would be evidence for QGP 
formation. In order to obtain the A/p ratio at the initial 
stage, one has to postulate from measurements at the 
final freeze-out stage including the nuclear absorption ef- 
fects. To this end, we use the Ultra-relativistic Quantum 
Molecular Dynamics (UrQMD) model [22| to simulate 
Au+Au collisions at ^s NN — 4.9 GeV and record the 
initial production abundances of antibaryons and high- 
mass antibaryon resonances as well as the final freeze-out 
abundances. We chose UrQMD because, as we discuss 
later, (1) it has been reasonably successful in describing 
many of the experimental results on hadron spectr a, as 
well as the average baryon density and the baryon emit- 
ting source size which are the essential ingredients for 
nuclear absorption, and (2) it does not have a QGP state 
or mechanisms mimicking a QGP state. Our strategy is 
then to compare the final freeze-out (A + S° + l.lT,~)/p 
ratio to data and use the initial A/p information from 
UrQMD to conjecture about what the data might be 
telling us. We shall demonstrate that the observed large 
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(A + E° + l.lS _ )/p is consistent with strong annihilation 
of A's and p's in the nuclear medium created in Au+Au 
collisions. Therefore, the data does not lead to a conclu- 
sion that QGP is formed. 

The UrQMD model has been applied successfully to 
explore heavy-ion reactions from AGS energies (Ei a b 
1 - 10 AGeV) up to the full CERN-SPS energy (E lab = 
200 AGeV). This includes detailed studies of ther- 
malization (23|, particle abundancies and spectra 0, 
HH, strangeness production (2(|, photonic and leptonic 
probe s l27l[28t , J/<3>'s [HI, and event-by-event fluctua- 
tions daliT 

UrQMD is a microscopic transport approach based on 
the covariant propagation of constituent quarks and di- 
quarks accompanied by mesonic and baryonic degrees of 
freedom. It simulates multiple interactions of ingoing 
and newly produced particles, the excitation and frag- 
mentation of color strings, and the formation and decay 
of hadronic resonances. The leading hadrons of the frag- 
menting strings contain the valence-quarks of the origi- 
nal excited hadron. In UrQMD they are allowed to in- 
teract even during their formation time, with a reduced 
cross-section defined by the additive quark model, thus 
accounting for the original valence quarks contained in 
that hadron [22J. 

Within the model, antibaryons are produced through 
string fragmentation. The Field— Feynman fragmenta- 
tion mechanism [32] , which allows the independent string 
decay from both ends of the string, is used in the UrQMD 
model [22]. The string break-up is treated iteratively: 
String — ¥ hadron + smaller string. The conservation 
laws are fulfilled. The essential part of this mechanism 
is the fragmentation function which yields the probabil- 
ity distribution p(z^ action ). This function regulates the 
fraction of energy and momentum given to the produced 
hadron in the stochastic fragmentation of the color string. 
For newly produced particles the Field-Feynman func- 
tion 13211: 
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is used. P(z) drops rapidly with increasing z. Therefore, 
the longitudinal momenta of the produced antibaryons 
are small; they are mostly produced in the central rapid- 
ity region with high baryon densities [33j . 

At AGS energies, antibaryon production is very rare. 
In order to increase statistics, we modified the string 
fragmentation routine in such a way that every string 
fragmenting process is repeated to a maximum of 1000 
times or until at least one antibaryon or antibaryon res- 
onance is produced. This drastically increased the ab- 
solute abundance of antibaryons, but does not alter the 
relative abundance among antibaryons. 

Once produced, an antibaryon may or may not annihi- 
late with baryons in the collision zone. The pp annihila- 
tion cross-section is well measured. Figure [1] depicts the 
annihilation cross-sections together with the elastic and 
total cross-sections [3J] . In this study, the pp annihilation 
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FIG. 1: (Color online) Parameterization (dashed curve) to 
the measured pp annihilation cross-section (open circles) 
along with those of elastic and total cross-sections. The 
parametrization is used in the UrQMD model. The Ap anni- 
hilation cross-section used in the model is deduced from that 
of pp by the Additive Quark Model. 



cross-section is given by 



T total / 



(2) 



where ^/s is in GeV. The total pp cross-section is taken 
from the CERN-HERA parametrization (34j (shown in 
Fig. HI . The other antibaryon-baryon annihilation cross- 
sections are, however, not well measured. UrQMD ap- 
plies a correction factor, given by the Additive Quark 
Model for these annihilation cross-sections 1331: 



1-0.4^1 fl-0.4^ 
3 A 3 
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Here sg is the strangeness number of the antibaryon and 
the baryon, respectively. There is a reduction of 13% due 
to each strange or anti-strange quark. For instance, the 
Ap annihilation cross-section is 



°Ap( Vs) = 0.87a Pp (yG) 



(4) 



Note that the relationship in Eq. (j4|) is for the same 
c.m.s. energy yfs of the Ap and pp systems. Generally, 
the Ap c.m.s. energy is larger than that of pp in heavy- 
ion collisions, so the reduction factor is lower than 0.87. 
The average c.m.s. energy square of a pair of particles in 
a chaotic system is approximately s w ra\ + m\ + 2E\E2, 
where rrii and Ei are the rest mass and total energy of 
the particles. The transverse distributions of p and A 
have been measured at rapidity about 1.2 [2lll3a|. Taken 
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the mid-rapidity transverse distributions to be similar, 
the average energies of the p and A at mid-rapidity are 
roughly 1.18 and 1.42 GeV, respectively. The c.m.s. ener- 
gies of pp and Ap pairs are therefore 2.13 and 2.34 GeV, 
respectively. This difference in y/s introduces an addi- 
tional reduction of the Ap annihilation cross-section by 
approximately 23% as opposed to the pp annihilation 
cross-section. Therefore, the effective relationship be- 
tween the pp and Ap annihilation cross-sections at the 
AGS is 



0.67(cTp P ) 



(5) 



Figure [5] gives an idea about the magnitude of the ab- 
sorption effect in Au+Au collisions by plotting the ratio 
of freeze-out A (or p) over that at the initial production 
stage. The ratio can be viewed as the "survival prob- 
ability" of A (or p) from initial production to the final 
freeze-out. The ratio is higher than one in forward and 
backward rapidities because the final dN/dy distributions 
can be broader than the initial ones. The left upper plot 
shows the "survival probability" as a function of rapid- 
ity, and the left lower plot shows that as a function of 
Pt for central Au+Au collisions. It is clearly seen that 
the largest absorption is in the mid-rapidity and low px 
region. The right plot shows the "survival probabilities" 
of mid-rapidity (\y\ < 0.4) and low px (px < 0.3 GeV/c) 
A's and p's in solid symbols and integrated over whole 
phase-space in open symbols, as a function of impact pa- 
rameter (b) in Au+Au collisions. As seen from the plot, 
about 99.9% and 99% of the mid-rapidity and low px p's 
and A's produced in central Au+Au collisions are anni- 
hilated. In other words, only 1 out of 1000 p's and 1 out 
of 100 A's in this kinematic region survive to freeze-out. 

On the other hand, if all the A's and p's are counted, 
then the "survival probabilities" are much higher, and 
are roughly constant over a wide range of impact param- 
eters in central collisions. However, this has implications 
on the interpretations of the measured absolute p yields. 
The measured yields (within a fixed rapidity window) has 
a less than linear increase with the total number of partic- 
ipants in Au+Au collisions [35[ • The power factor of the 
increase is 0.74. If we take into account the absorption 
effect shown in the open symbols in Fig.^c), then the re- 
stored initially produced p's would have a stronger than 
linear increase with the total number of participants. The 
power factor would be about 1.46. 

Now back to Fig. [2] In the simple picture of a sphere of 
baryons with a uniform density p and radius R, the sur- 
vival probability of an antibaryon produced at the center 
is 



-Psurv = cxp(- a pL) , 



(6) 



where a is the annihilation cross-section. By using the 
factor 0.67 from Eq. ([5]) as our power factor, we have 

(7) 



p_ __ p0.67 

A.surv p,surv ' 



and obtain the dashed curves in Fig. [5] given the p data 
points. We find good agreement with the calculated A 
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FIG. 2: The ratio of number of A's (open symbols) and p's 
(solid symbols) at final freeze-out over that at the initial string 
fragmentation stage, calculated by UrQMD for central (b < 1 
fm) Au+Au collisions at y/s NN = 4.9 GeV, as a function of 
(a) rapidity but integrated over px, and (b) px but integrated 
over rapidity, and for all (minimum bias) collisions, as a func- 
tion of (c) impact parameter b for the mid-rapidity (\y\ < 0.4) 
and low pr (pr < 0.3 GeV/c) region (circles) as well as for in- 
tegrated over whole phase-space (squares). At AGS energies, 
antibaryons can be produced only at the very early time from 
string fragmentation, and are then annihilated by baryons at 
later times. However, the ratio can be larger than one be- 
cause the final freeze-out rapidity and/or pr distribution can 
be broader than the initial distributions. With this caution, 
the ratio may be viewed as the "survival probability" of A and 
p through the nuclear matter created in the collisions. The 
dashed curves indicate a simple optical model prediction of 
the A survival probability using the p result and the different 
annihilation cross-sections of A and p with nucleons. See the 
text. 



survival probability except for the forward and backward 
rapidity regions. 

Figure [3] shows the freeze-out ratio of (A + E° + 
l.l£ - )/p at mid-rapidity (|y| < 0.4) and low px (pt < 
0.3 GeV/c) in Au+Au collisions at AGS energy as a 
function of the collision impact parameter b. The ex- 
perimentally deduced ratio is reproduced from Ref. (l6j 
with b values obtained from the centrality bins. The 
UrQMD ratio is in a good agreement with the data. 
The (A + S° + l.lfr )/p ratio of the total particle yields 
are also shown. It is clear that the large ratio in the 
mid-rapidity and low px region in central collisions is 
largely due to the kinematic cut. Note that the total 
yield (A + E° + l.lE~)/p ratio in p + p interactions as 
calculated by UrQMD is in a good agreement with the 
trend from the heavy-ion results. 

As UrQMD successfully describes the data, it is inter- 
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FIG. 3: UrQMD calculation of the freeze-out (A + E° + 
l.lS~)/p ratio at mid-rapidity (|y| < 0.4) and low pr (pt < 
0.3 GeV/c) as a function of impact parameter b (solid tri- 
angles) in Au+Au collisions at ^/s NN — 4.9 GeV, compared 
to experimental results in similar kinematic region (solid cir- 
cles). Note that UrQMD can reasonably reproduce the data. 
The UrQMD total yield (A + E° + l.lE")/p ratio (no kine- 
matic cut) is shown as open triangles. The UrQMD total yield 
(A + E° + l.lE~)/p ratio in p + p collisions at = 4.9 GeV 
is indicated by the box drawn at b = 14 fm. 



esting to examine the ratio at the initial production stage. 
Since high-mass antibaryon resonances are present at the 
initial stage, it is important to correct for the feed-down 
contributions of these resonances to p's, A's, and E's. We 
find that the initial ratios of (A + S° + l.l2 _ )/p « 0.46 
and A/p ~ 0.14 are below unity and are roughly indepen- 
dent of centrality. The initial A/p ratio is consistent with 
that calculated in isospin averaged nucleon-nucleon colli- 
sions. The centrality independence can be readily under- 
stood because the string fragmentation function does not 
know about the centrality of the collision and there are 
no other antibaryon production mechanisms in UrQMD 
that are centrality-dependent (e.g. due to QGP produc- 
tion) . 

As UrQMD well reproduces the freeze-out values for 
the (A + S° + l.lX~)/p ratio, it is reasonable to suggest 
that the data indicates no large ratio of A/p at the initial 
stage, and therefore does not imply QGP formation. 

The (A + fl )/p ratio has been also measured in heavy- 
ion collisions at the SPS. The data show a decrease of the 
ratio from central to peripheral colllisions (36|. The de- 
crease was also observed in UrQMD where the A/p ratio 
drops rapidly with increasing b from 1.3 to 0.5 |37l |. This 
suggests an interplay between particle production and 



FIG. 4: The excitation function of the (A + E°)/p ratio 
from UrQMD in the AGS, SPS, and RHIC energy range. 
The UrQMD calculations were done for \y\ < 0.4 and px < 
0.3 GeV/c in central Au+Au collisions (6 < 1 fm). 



subsequent annihilation also at the SPS energies. In pe- 
ripheral (large b) collisions the A production is basically 
the same as in p + p reactions. A's and p's are produced 
via the fragmentation of color flux tubes (strings). The 
production of (anti-) strange quarks in the color field is 
suppressed due to the mass difference between strange 
and up and down quarks. This results in a suppression 
of A over p by a factor of 2 (A/p w 0.3-0.5 in p + p). 

The (A + Ti°)/p ratio in high energy heavy-ion colli- 
sions at RHIC was measured to be singinificantly smaller 
than that at the AGS and SPS energies [H, H| . This is 
an experimental demonstration of the importance of the 
net-baryon density on the (A + S°)/j> ratio. At the top 
RHIC energy, the net-baryon density is small at mid- 
rapidity [40l . |41| . resulting in an insignificant effect on 
the (A + S°)/p ratio. The (A + S°)/p ratio measured 
at RHIC, therefore, reflects more truly the initial pro- 
duction stage. This also demonstrates that a mere large 
(A + ratio is by no means a signature of the QGP 

formation, as the QGP formation is much more likely at 
RHIC than at the AGS. One has to take into account the 
effect of antibaryon absorption before connecting a large 
(A + TP)/p ratio to the QGP formation. 

In Fig. H we show the (A + t°)/p ratio from UrQMD 
as a function of the collision energy y/s NN . Here we have 
kept the same parameters as for the previous AGS anal- 
ysis for consistency, i.e. a cut on mid-rapidity and low pt 
and the same centrality. The ratio is found to decrease 
with y/s NN because of the decreasing net-baryon density 
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with increasing collision energy. 

In conclusion, a large A/p ratio is observed in central 
Au+Au collisions at the AGS. A strong increase of the 
A/p ratio from peripheral to central collisions is indi- 
cated from the experimental data. The hadronic trans- 
port model, UrQMD, can satisfactorily describe the large 
ratio in central collisions and the centrality dependence 
of the ratio. According to the model, the experimentally 
deduced large ratio of (A + S° + in the mid- 

rapidity and low px region is mainly due to the strong 
and different absorption of these antibaryons. The mea- 
sured large A/p ratio in itself is by no means a signature 
of the QGP formation. An increase of the A/p ratio from 
peripheral to central collisions was also observed at the 
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